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Abstract. We study the asymptotic expansion of the smeared L'^-trace of 
/ e~'^ where P is an operator of Dirac type, / is an auxiUary smooth smear- 
ing function which is used to localize the problem, and chiral bag boundary 
conditions are imposed. Special case calculations, functorial methods and the 
theory of and r;-invariants are used to obtain the boundary part of the 
heat-kernel coefficients ai and 02. 



1. Introduction 

Local boundary conditions for operators of Dirac type have been studied in the 
physical and mathematical literature with a variety of motivations over many years. 
Some key points in this respect are as follows. 

(i) Local boundary conditions for massless fermionic fields ruled by a Dirac operator 
can be applied to one-loop quantum cosmology |12| . |18j and are part of the inves- 
tigation of conformal anomalies {221 in Euclidean field theory jJH!- Moreover, they 
are the first step towards analyzing boundary counterterms in supergravity theo- 
ries, with the associated unsettled issue of proving finiteness ^21 or lack of finitcness 
|14| of supergravity theories on manifolds with boundary. In other words, the local 
boundary conditions for fermionic fields arc part of a general scheme |28) leading to 
locally supcrsymmctric boundary conditions for fermionic and bosonic fields |12j . 
P] , and hence can be used to test perturbative consistency of supergravity models 
in cosmological ^Hl, 1201 or field-theoretical backgrounds. 

(ii) Local boundary conditions of chiral bag type arc a substitute for introducing 
small quark masses to drive the breaking of chiral symmetry j3(J) . One of the first 
papers where the chiral boundary conditions were introduced is the work by Hrasko 
and Balog [221 j and one of the first applications to chiral bag models is presented 
in 122. 

(iii) Chiral bag boundary conditions have been recently proved to lead to a strongly 
elliptic boundary- value problem for the squared Dirac operator jS], and the associ- 
ated global heat-kernel asymptotics has been investigated in detail, on the Euclidean 
ball, in [21]. An early paper on the role of boundary conditions for Dirac operators 
is in the framework of fermionic billiards [2], studied even earlier by Berry and 
Mondragon [Jj. 

For more general Ricmannian manifolds with boundary, however, the investiga- 
tion of such a global asymptotics in the chiral bag case is, to our knowledge, an 
open research field, and it appears desirable to understand how far can one go by 
exploiting functorial methods (e.g. conformal rescalings of the metric) and special 
case calculations, which are tools frequently used in invariance theory |23| . |27| . 
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For this purpose, both algorithms are exploited in our paper, where the general 
mathematical setting is as follows. 

Let m = 2m be even and let P = 7jVj + ■0 be an operator of Dirac type on a 
compact oriented Riemannian manifold M of dimension m, where V is a compatible 
unitary connection, i.e. V7 = 0. The spinor space has then dimension ds = 2™, and 
the 7-matrices can be taken to be skew-adjoint and obeying the Clifford relation 

Near the boundary, let em be the inward-pointing unit normal and 7^ be the 
projection of the 7-matrices on e^. Moreover, the generalization of 75 to arbitrary 
even dimension is provided by 

(l.a) 7 = f"7i...7„. 

The squared Dirac operator is studied with local boundary conditions of chiral bag 
type. These boundary conditions involve a real angle 6 and they read 

(l.b) ^-^\aM - 0, 

where we have introduced the 'projectors' 

(1-c) ^ i (1 ± e''^77m) . 

Under the above assumptions, the squared operator is an operator P of Laplace 
type [231 ■ The associated heat kernel can be defined as the solution, for t > 0, of 
the heat equation 

(l.d) (^^^+P^U{x,x';t) = Q, 

obeying the initial condition 

(l.e) lim / dx'U{x,x;t)'^[x) = '^{x), 

jointly with the boundary conditions Be defined by 

n_C/(x,rr';0L6aA/ = 0^ n_F,[/(x, x'; i)LeaM = 0. 
Here, dx' denotes the Riemannian volume element of the manifold M and P^ denotes 
the Dirac operator with respect to the variable x. The L^-trace of the heat semi- 
group is obtained by integrating the fibre trace Try of the heat kernel diagonal 
U{x, x; t) over M, and reads as 

(l.f) Ti-L2(e~*^) ^ I dxTrvU{x,x;t). 

In our paper, following we are interested in a slight generalization of the 
previous equation, where e~*^ is 'weighted' with a smooth scalar function / on M. 
More precisely, we are interested in the asymptotic expansion as i — > 0^ of the 
functional trace 

(l.g) Tri2 (/e-*^) = dx f{x)TTvU{x, x; t). 

The results for the original problem are eventually recovered by setting / = 1, but 
it is crucial to keep / arbitrary throughout the whole set of calculations, as will be 
clear from the following sections. 

The asymptotic expansion of such a functional trace has the form 

00 

(l.h) Tr^. (/e-*^) ^ P, Be). 

ri=0 

Note that there is a change of convention in the indexing of the Seeley coefhcients 
with respect to the work in [23, i-S- our a„ is denoted therein by a„/2- The 
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coefHcients an{f,P,Be) consist of two different parts, the interior part a^'{f,P) 
and the boundary part a^^^(/, P, Bg), i.e. 

(1 -i) aM,P,Be) = a^' {f,P) + af ^ {f,P,Be). 

The interior parts a*^ (/, P) are obtained by integrating some geometric invariants 
(see below) over M and are independent of the boundary conditions. By contrast, 
the boundary parts a'^'^^ {f, P,Bg) are obtained by integrating some geometric in- 
variants over the boundary dM and these parts depend in a crucial way on the 
boundary conditions. They will be the main concern of our research from now 
on. The interior invariants are built universally and polynomially from the metric 
tensor, its inverse, the Riemann curvature of M , the bundle curvature (if a vector 
bundle over M is given) and the endomorphism (or 'potential' term) in the squared 
operator P^. By virtue of Weyl's work on the invariants of the orthogonal group, 
these polynomials can be formed by using only tensor products and contraction of 
tensor arguments. Here, the structure group is 0(m), m being the dimension of M. 
However, when a boundary occurs, the boundary structure group is 0{m — 1), and 
the Weyl theorem is used again to construct invariants. 

The structure of the article is as follows. In Section 2 we write down the general 
form of the leading coefficients oi and a2. The special case calculation of 
and different functorial techniques are used to determine part of the numerical 
multipliers of the geometric invariants. Further special cases are shown in Section 
3 and the complete ai and 02 coefficients are determined. We end the paper with 
concluding remarks. 



2. Determination of the leading coefficients 

We first write down the general form of the leading two boundary contributions to 
the heat kernel (hereafter, Laa is our notation for the trace of the extrinsic-curvature 
tensor of the boundary). 

Lemma 2.1. Let f be scalar. There exist universal constants Ci{9,m) such that 
(hereafter our notation for the invariant integration measure on dM is simply dy ) 

(2.a) ar(/,P,6e) = (4^)-('"-i)/2 / dy Try(ci(0, m)/). 



dM 



af^(/,P,60) = (4^)-'"/2 / dyTrv{c2i0,m)Laaf + C3{0,m)f^jj„ 



dM 

(2.b) +€4(0, m^fil^^ryi + 05(0, m)fi'j + ce{9, m)fil: + 07(6*, m)f.,n) ■ 

Proof. This is a direct consequence of the Weyl theorem on the invariants of the 
orthogonal group [23, as we said at the end of Section 1. □ 

We next determine the universal multipliers Ci{9,m), i = 1,...,7. We first ex- 
ploit a known special case. As usual, the hypergeometric function is denoted by 
2Fi(a,6; c;z). 



Lemma 2.2. We have 



(2.c) ci {9, ™) = ^ (cosh™-! 9-1) 



4 



GIAMPIERO ESPOSITO, PETER GILKEY AND KLAUS KIRSTEN 



Proof. In [2] the heat kernel coefRcients for the given setting have been evaluated 
on the Euclidean ball for the case ip — and / = 1. The results obtained were 

a, ^ ^pl^ (cosh™-^ - 1) , 
2™r(m/2) ^ ^ ' 

(2m — 5)ds ds 

a2 - 



3-2'"r(m/2) 2"T(m/2) 
iFi (^l,^^;i;tanh2^?^ - (m - 1) ^Fi (^1, | tanh^ 

The volume of the sphere, which is the boundary of the ball, is 

Using this to rewrite the coefficient ai shows the assertion 1)2. c|l . which agrees with 
Eq. (40) in for the first boundary correction to the partition function. 

To show (|2.d|l we first use the Gauss recursion formula, see e.g. |2S|, equation 
9.137.17, 

7 2F1 (a, /?; 7; z) - (7 - /3) 2F1 (a, /3; 7 + 1 ; z) - /? 2i^i (a, /? + 1; 7 + 1 ; ^) = 0, 

to write 02 for the m-ball as 

ds f 2m — 5 ,„ , ^ m — 1 3 , y „ 

02 = 7^-^ < h 2 - m) 1, ; -;tanh^6' 

2"r(TO/2) \ 3 ^ 2'2' 

Comparison with the general form (|2.bp then shows assertion Ij2.dp . Note that in 
the given setting, i.e. with = and / = 1, the C2{0,m)Laaf term is the only 
term contributing. □ 

Remark 2.3. For 6 = the boundary conditions reduce to standard boundary 
conditions of mixed type. For 9 = we have 

ci(0,m) = 0, 

1 (2m-5 , ,1 1 1 - TO 1 

C2 0,TO = ^^^^+2-m.l =- —^ = - 

2(m — 1) t 3 J 2(to— 1) 3 6 

To achieve comparison with the known results for mixed boundary conditions note 
that the auxiliary Hermitian endomorphism x needed to define the splitting of the 
spinor bundle is |27j 

X = -11m- 

Let 

n± = \{i±x) 

be the projection on the ± eigenspaces of x- Mixed boundary conditions are then 
defined as 

Bip = n_vj © (v,„ + S')n+(/9|g^^ = 0. 

The relevant S for the given setting is 

S = ——LaaTl+. 

Using the fact that Ttv{x) = 0; Tr\/(n-|-) = ds/2, the coefficients for the relevant 
mixed boundary conditions follow to be 

af'\l,P,Bo) - 0, 

af^(l,P,So) = (4^)-"/"^ I Try(2L,a + 125)d2/=(4^)-"/2i J Trv{-Laa)dy, 

dM dM 

in agreement with our findings for ci(0,to) and 02(0, to). 
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We next exploit the fact that the connection V is not canonicaUy defined. To 
simphfy the notation shghtly we assume a locahzing function / = 1. 

Lemma 2.4. We have 

Ci{9, m) = 0. 

Proof. Let a.^ be a skew-adjoint endomorphism of the spinor bundle commuting 
with the Clifford structure 7, [CTi,7j] = 0. Then 

defines a smooth one-parameter family of compatible unitary connections. We define 
to ensure that 

P(e) = 7«V.(e)-fV'(e) = ^ 
is unaffected by the perturbation; the boundary condition also remains unchanged. 
Therefore, the heat trace coefficient (j2.b(> remains unchanged. Using 77^ = ~7i7 
we evaluate the variation S = {d/de)\^^o of the single terms for (Tq = 0, am 7^ 0: 

STTv{c2{0,m)Laa) = 0, 

(5Try(c3(6',m)V'77m) = - Tr\/(c3(6', m)7„,(T„77,„) = Try (03(6*, 771)0-^7) 
= -Try(c3(6', m)CT™7„77™) = - Try (03(6', m)(T,„7) = 0, 

(5Try(c4(6',TO)V'7m) = - Try (04(6', m)7™(T„7„J = Try (04(6', m)cr™), 

(5 Try (05(6*, to) V'i') = -Try(c5(6', m)7™cr™7) - Try (05(6*, m)cr™7) 
= Try(c5(6',TO)cr„7™7) = Try(c5(6', TO)7,„cr„i7) = 0, 

6TTvice{9,m)ip) = - Try(c6(6', m)7™(T„J = - Try (c6(6', m)7™cr„77) 

^ Try(c6(6',TO)77™crm7) = Try (c6(6', m)7™cr„i) = 0. 

For the coefficient to remain unchanged we need 04(6, m) ~ 0. 

Considering ctq 7^ and cr„, = does not produce any new information. □ 

To find more information about the remaining unknown multipliers, one might 
enlarge the setting and allow for an endomorphism-valued /. However, apart from 
the fact that the number of invariants goes up to 36 and the calculation gets cum- 
bersome, this does not produce any relevant information for our problem and we 
do not present further details. 

Instead, we next exploit conformal rescaling techniques. 

Lemma 2.5. We have 

1 



07(6', m) = - 02(6*, TO, , 

TO — Z \ O 

Proof. Let / be a smooth function with /|aj\f = 0. Define ds'^{e) := e^'^^ ds^ and 
P(e) := e~'^^P. Let V be a compatible unitary connection. We expand P = 
^'^'S/ + tjj with respect to a local coordinate system x = (xi, ...,Xjn) and use the 
metric to lower indices and define 7^. If we define 

V(e)a, Va^. + ^e(./;.7''7M + U)^ 
results of |^ show that V(e) is a compatible unitary connection. Furthermore, 

Note that the boundary condition remains unchanged under conformal variation. 
The heat kernel coefficients satisfy the equation 

^1, Pie), Se) = (to - n)a„ (/, P, Be 



e=0 
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To study the numerical multiplier c-j{9, m) we need the variations 
d 



de 

d 



T(e) = -2/t - 2(771 - 1)A/, 

e=0 



de 



Laaie) = -fLaa - (m - !)/;,„, 

where r = Rijji is the scalar curvature. Applying equation (j2.e|l shows the assertion. 

□ 



Remark 2.6. Note that, despite the appearance, the multiplier cj{9,m) is well 
defined in dimension m = 2. Using the result for C2{0,m) given in equation 
we obtain explicitly 

(2.f) cj{e,m) = -Ul- 2^^! fl,!^;|;tanh2( 



2 [ - - \^ ' 2 '2' 
For 9 — this agrees with the previous computation for mixed boundary conditions. 

3. Relating the zeta and eta invariants 

In order to determine the numerical multipliers 03(0, m), 05(0, m) and ce{9, m) we 
relate the zeta invariant to the eta invariant. We will then evaluate the eta invariant 
on the TO-dimensional cylinder and ball for the case of an endomorphism- valued /. 
On the ball we will restrict to the choices / = 1 and / = 7, respectively, which 
will allow us to find 05(^,771) and C6(^,77t,). Instead, on the cylinder we can deal 
with general /. Performing the two special case calculations is strictly speaking not 
necessary, but provides helpful crosschecks of the answers obtained. 

To distinguish the coefhcients in the heat trace, Tri2(/e~*^'), and in the trace 
related to the eta invariant, Tr]^2[fPe~*^ ), in this section we use the notation 



n 

Tr^. (fPe-'P^-) ~ if, Be) ■ 



The result we are going to need is the following: 

Lemma 3.1. Let f e C°^{End{V)) and let P, := P + ef. We then have 

d,ai{l,PlBg) = -2al_,if,P,,Bg). 

Proof. The proof is insensitive to the boundary conditions imposed and parallels 
the proof in [ini. □ 

Remark 3.2. The very useful property of this result is that the coefficient for 
the zeta invariant is related to the coefficient a^ -^ for the eta invariant, which will 
have a significantly simpler structure. 

In order to apply Lemma 13.11 to the coefficient 03 we need the general form of 
the coefficient. 

Lemma 3.3. Let f G C'°^ {End{V)) . There exist universal constants di{6,m) such 
that 

a^,'''\f,P,Be) = 

(4^)-m/2y Tw {di{9,m)f + d2{9,m)f^ + d3{e,m)f-frn + di{e,m)f^-fm} . 

dM 

Proof. This follows immediately from the theory of invariants taking into account 
that / is in general a matrix-valued endomorphism. □ 
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Remark 3.4. Lemma \c/.1\ relates the universal constant dj{9,m), j — 1, ...,4, with 
Ci{9,m), i = 3, ...,6. In particular we have 

C3{9,m) = — 2^4(6*, m), €4(0,771) = — 2^3(6*, m), 
c^{e,m) ^ ~2d2{9,m), CQ{e,m) = -2di{e,m). 

From Lemma ^2.4\ we conclude d'i{6^m) ~ 0. We evaluate di(9,m) and d2{0,m) 
for the example of the ball and thus find c^{9,m) and CQ{6,m). We also evaluate 
di{6, m), d2{0, m) and £^4(0, m) for the example of the cylinder. This provides checks 
of the answers for c^{6,m) and CQ{9,m) and in addition determines di{9,m) and 
thus C'i{6, m). 

For the ease / = 1 we proeeed as described in The case / = 7 is based 

upon this calculation and therefore we need to present some details for the case 
/ = 1. We first summarize properties of the spectral resolution for the Dirac 
operator on the ball. Let P = 7iVj; be the Dirac operator on the ball and let us 
denote by (p± its eigenfunctions obeying the eigenvalue equation Pip± — ±^ip±. 
On writing the eigenvalue equation in this form we have fj, > 0. Later on we will 
write the eigenvalues of P as A = ±/i, such that |A| = Modulo a suitable radial 
normalizing constant C, we may express |15) 

Here, Jy{z) are the Bessel functions and Zj_"''(ri) are the eigenspinors of the Dirac 
operator P on the sphere 

Pzi"'(r!) =±(^n+ Z^^\n) for n = 0, 1, ... 

The degeneracy dn{m) for each eigenvalue is 

dn{m) dimZ^"''(f2) = ^4 ^ ^ ^ " ^ 

We next apply the boundary operator which reads explicitly, from Eq. (l.b), 

1 / 1 -le" 

2 V ie-'^ 1 

to the solutions H3.a|) and IjS.bp . This produces the following eigenvalue conditions: 

(3.c) J„+™(A*)TeV„+™_i(/^) = for 

(3.d) J„+^(/i)±e-«J„+^_i(/^) for ip[^K 

These equations allow us to rewrite the eta function 
77(s;l,P,6e)=^sgn(A)|A|-^ 

A 

in terms of a contour integral and to apply the techniques described in detail in 
ini El m ■ The coefficients in the asymptotic expansion for the eta invariant are 
then determined by evaluating residues of rj according to |23 | 

(3.e) Res r^{m n; 1, P, Bo) - ^""/^-rV^f ■ 
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For notational convenience wc introduce p = n + m/2 — 1. Starting point of the 
analysis is [22 



1 /* r1 i 

t'o 27rz7 dk 1 



Jp+i(fc) 

^ Jp{k) 



where F is a suitable counterclockwise contour enclosing all solutions of the equa- 
tions H3.c|l and Ij3.d(l . After deforming the contour to the imaginary axis this gives 



TTl \ 2 J J dz T I jpO£p±l 

n=0 Q ± T IL r ^ 



(£) 

z /p(z) 



where in the last step we have used the recursion for the modified Bessel function 

m 

(3.f) /,+i(z) = /:(z)-^/,(z). 

In order to recover the coefficient al we only need to consider the leading term in 
the uniform p —^ oo asymptotic expansion of the Bessel function , 



/;(M (i + fc^)v^ (i\c(^ 

Ip{kp) k \ \p 



Hence we only need to find the residue of 



Ms; 1) = - cos (f ) ^ d^p- y dkk--ln ^ - iO ^ ^0) 

at s = m — 1 . 

We first observe that the summation over n produces a multiple of the Barnes 
zeta function which is defined by 

In detail we have 

^ dn{m)p^'' -dsCs (s, y - 1 



In order to perform the fc-integral we first combine 9 and —9 and evaluate the 
logarithmic derivative to give 

d 1 + ^ (1 - VI + fc2) 4isinh6l 



dfc + VTTA?2(2 + /c2 + fc2cosh(26i)) 
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The relevant fc-integral therefore reads as 

CO 

dkk-'{i + k^)-^''^— — p — 5- 



(cosh(20) - 1) 







From here, with |25j 

'1 - s 



r 



cos(f)r(i±^: 



we easily compute 



m 



^o(s;l) = -^^Y^lliysinheai^i (^l,^;2;-sinh2 0jCB(s,y-l 
On using [23 

r(m-i) 2'"-2r(f)' 

the coefficient can be cast into the form 

<(l,P,Se) = ir (y) Res 77(771 -l;l,P,Se) = ir(y) Res ^o(m-l;l) 
= - sinh ^ ~ ^, , 2F1 (1, 1 - ^; ^; - sinh^ 6> 



Comparing this with the answer on the ball expected from Lemma 13.31 
a1{l,P,Bg) ^ (4^)-™/2vol(5™-i)4di(^,™) = ^„^^^^^^ rf,di(g,m), 



™~1 TTlS ,9 



we read of 

(3.g) = -^^^sinh0 2P1 (^1,1- y;§;-sinh 

From Remark 13.41 we then get 

cq{9, m) = (m — 1) sinh 6 2F1 ^ — ^ sinh^ ( 

To find the universal constant 05(9, m) we perform the calculation on the ball with 
/ = 7. This choice complicates the analysis significantly because the normalization 
constant C and further integrals over products of Bessel functions come into the 
play. First we note that if ri{s]x,y) denotes the local eta function, then 

77(5; x,y)^Yl ^ (^)* ^+ ^ (y^ - (^)* ^-^ (y^} ■ 

We want to analyze 

(3.h) Tr (777(5; = E/^"' {< > - < ¥^^^^17^^^^ >} , 
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with < (pi\ip2 > denoting the Hilbcrt space product 

< fl\f2 >= / dx Lpl{x)ip2{x). 



M 

Since 

■ 1 



7 = 



-1 



changes the sign of the lower chirahty, the normaUzation constant C docs not cancel 
in the Hilbert space products appearing in l|3.h|) . but instead values of various 
integrals occur explicitly. We first observe that 



1 

1 



We use 



= dr r{Jp^^{iir) + Jp{^ir)). 



dr rJl{^r) = i - Ji,_i(/x) Ji,+i(/i)} 



to find 



We use the implicit eigenvalue equations (|3.c() and (|3.d|l together with recursion 
relations for the Bessel functions 

Jp+2ip-) = ^^^^-^ Jp+i(/i) - Jp(/i), Jp-l{^-l■) = —Jpip) - Jp+i{^), 

to simplify the normalization constants Cj_^'' for the different spinors iy9±^''. We 
obtain 

^(±) ^ I 

(^(±) ^ I 

Proceeding in the same way for the quantities < 'p'^^^ip^^^ >, we find 



e±« f f 



/i + /ie±2f zp(2„+l)e±e 2 cosh 61 ,, zr £±1/2 



2 cosh 6* -u £±1/2 ■ 

^ cosh e 

Using these results in ()3.h|l . we obtain the following contour integral representation: 

n=0 p coshe 

(3.i) 

^ ' 4-Kicoshe^ ^ V I P+i/2 



Note that the counterclockwise contour must only include the zeroes of the equations 
(|3.c|) and (|3.d|l such that the appropriate summation over eigenvalues results. The 
poles at k = {p + 1/2)/ coshfl should lie outside the contour because they have 
been introduced by the normalization integral and need not be summed over. The 
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situation is similar to the analysis for radial smearing functions, see jl7) for more 
details. This observation is important because when shifting the contour towards 
the imaginary axis additional contributions result. Using the index p for all Besscl 
functions an intermediate result reads as 



1 

27ri cosh 6 
1 



■ cos 



n=0 



2TTi cosh 
1 



cos C^j ^ dn{m) / dz 



n=0 



-s d 



■ , p+1/2 
~^ coshe 

^In [/;(z)-f/p(z) + ze%(z)] 



2 cosh ( 



n=Q 



cosh 9 



^ 1 



cosh 9 



_£+l/2 



(3.j) 



The last contribution resulting from the shifting of the contour can be given in 
closed form by using the differential equation for the Bessel function |25| . 



We calculate 



-ln(j;(fc) + (e«-^)j,(fc) 



dz^ 



,e P 



Id ^ ^ 
z dz 



Mz) = 0. 



P+l/2 
" coBh e 



4'(fc) + ^J,(fc) + (e«-f) j;(fc) 



J^{k) + {eo - f) J,(fc) 



.._ p + l/2 



sinh 61 



„_p+l/2 



cosh0 
2p+l' 



Adding the contributions from 9 and —9 the sinh0 terms cancel and the summation 
over n leads to C,b{s + 1, (m — l)/2), which has no pole at s = m — 1. Therefore, for 
the present purpose this term is irrelevant. 

In the remaining integrals in we need, as before, only the leading term in 
the Debye asymptotic expansion of Bessel functions. Explicitly, with x = 1/ cosh 0, 
we obtain to leading order 

cos °° 7 



n=0 



k — ix k + ix 



2cos(f) 



TT cosh 9 



^dnim)p " / dk 



n=0 



fc2 + x2 



The /c-integral is 



dk 



k-Wi + k^ 1 ri 



fc2 



2x2 



1 - s 1 



and hence 
^o(s;7) = 



1 cosh6i r f 

2 ■ 



2i^l 1 



1 - s 1 



; -; - sinh 9 ] (b (s, — - 1 
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The residue is easily evaluated and via H3.c|l we compare it with the form given in 
Lemma [3.31 to read off 

(3.k) d2{e,m) = -^coshe 2F1 M,l-y;^;-sinh'( 



which implies 



05(6*, m) = cosh 6* 2F1 (1,1- y;i;-sinh^( 



Remark 3.5. In the calculation just described it is the argument — sinh^ 9 that 
occurs naturally in the hypergeometric functions. Instead, the constant €2(0, m) 
in i2.d!^) and 07(0, m) in have been given using tanh 9. In order to provide 

answers in a unified way one might use the transformation formula |25j 



2Fi(a,/3;7;z) = (1 - z)"" 2F1 I a,j~P;r, 



1 

to write 

C2{e,m) = 2(m-l) {^^^^^ + (2-"^)cosh'g 2^1 (^1, 2 - ^; ^; - sinh^ I 



C7(6l,m) = --|l-cosh'6l2i^i (^l,2-y;-;-sinh 

Remark 3.6. Note that, despite the complicated appearance of the universal con- 
stants, for each specific dimension m a simple function of m and 9 results. In 
particular, whenever the second argument of 2F1 is or a negative integer, the 
hypergeometric function reduces to a finite polynomial in sinh^ 9. 

In order to find the missing multiplier C3{9,m) we present a calculation on the 
cylinder. In order to summarize previous results [S] we need to provide some nota- 
tion. Let M = IR+ X be an even dimensional cylinder equipped with the metric 
ds^ ~ dx^y^ + ds^, where Xm is the coordinate in IR+ and plays the role of the nor- 
mal coordinate, and dsjq is the metric of the closed boundary N . The coordinates 
on N are denoted by ?/ = (yi,?/2, ■•■,ym-i)- To write down the heat kernel on M 
for = (7iVi)^ with boundary condition Bg, we call (j)uj{y) the eigenspinors of the 
operator B = 77m7aVa, corresponding to the eigenvalue uj, normalized so that 

with 6™^^ the Dirac delta function, and 

dy (PlivMy) 1. 



TV 

Finally we need a; = (y, a;,„), ^ = x„i-x'^, = Xm+x'^, Uu,{ri,t) = ^-A/iw tanh 6*, 
and the complementary error function 



2 [ 

erfc(a;) = —= I d^e 

We then have [S] 

U{x,x';t) = -^Y.cl,*{y')My)e-'^^*{(e"^-e^)l 

(3.1) +^^T^7^ + yM^tanh0e"'('''*)erfc(M^(77,i)) 

cosh" (6^) L 

(Note that although the formal appearance of the heat kernel is identical to the one 
in [S], equation (5.1), the meaning of 11+ is slightly different. The reason is that 



— n 

e 4* 
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[S] considers the boundary condition resulting from, in our present notation, n_j_ 
whereas we consider the one resuhing from n_. FormaUy the transition is obtained 
by reversing the sign of the normal and by using our present notation for II-i-.) As 
remarked in [S] the first term is the heat kernel on the manifold IR x iV, which does 
not encode any information about the boundary contribution. In the following, 
without changing the notation, we will ignore this term and we will determine the 
boundary contributions to the eta invariant from the remaining terms. 

Let / G C°°(End(y)), then we want to consider TTi^2(f[PxU{x,x';t)]x=x'); note 
that the derivatives need to be performed before the coincidence limit x = x' is 
taken. Given we have the local form of the heat kernel we can in principle deal with 
an arbitrary /. For our present purpose it is easiest to assume / = /(y) only such 
that the Xm-integration can be done without complication. 

It is natural to introduce the heat kernel UB{y,y'',t) of the operator B^, 

c/B(y,2/';t) = ^C(y')0-(yK"'*; 



furthermore, to make the single steps easier to follow we use the splitting 



m{x,x';t) = --==^</>*(y)</>^(y)e-"'*e^. 



U2{x,x']t) 



2n+n* 



'47rt ^ cosli (6*) 

1 + v/(7rt)t^tanh6le"'('''*'erfc(uc^(?7,i))j 



Acting with P and performing the ^m-integration, intermediate results are 



dXrfYi f^PxU\{tl^y ^ Xrn: X^-^^t^\xjn—x' — 



(3.m) 



■^=7;flmUBiy,y';t) - ]-fj,n^ByUB{y,y';t), 



dxm f[PxU2{y, 



1 



2 eosh^ e 
(3.n) 



f'ymUB{y,y';t)Il+Il* 



> tanh 9e 



1 



2 cosh' 

Here, we have used the relation 
1 d 



fjmWBiy, y'; t)U+n*^e 



* tu tanh 



^erfc(— Vttj tanh 9) 
'erfc(— -v/to; tanh 6*). 



2 dxr. 



xf^/t+U-^{2Xm,t) 



CTic{u^(2x„i,t)) 



— +tjtanh6'e"-(2^'-*)erfc(u„(2x™,t)) . 
irt 

Whereas the asymptotic t ^ behaviour in (|3.m() could be easily found from the 
corresponding (known) behaviour of the trace of Ub, the same is not as simple for 
the result in H^.nfl . We have found it most convenient to perform the L^{N)-trace 
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and to relate the above equations to the zeta and eta function via 

oo 

1 



r(.) 



dt f-'TTL2 fe 



tP- 



oo 

and to evaluate the asymptotic t —* expansion from (|3.e|) and 



(3.0) 



r(z) 



For H3.m|l the associated relation is readily found, 



mis-J,P,Bg) 



1 1 ri 



/4^2r(^)^ V2 
In order to proceed with H3.n|) we note first that 

tuj"^ tanh^ 6 - tuj"^ 



cosh^ e 



erfc(Vta; tanh 9) — 1 + erf(v^u; tanh 9), 
with the error function 



erf(.T) 

The resulting t-integral then is 



dt e"* 



dtt'^^e colh2 9 fi + cTf{VtLL> tanh I 



cosh'* 



CO 



s+1 



r 



s + l 



-^T 1 + - sinhfegn(w) 2Fi(-,l + -;-;- sinh 



1 



s 3 



2 2' 



This produces the following contributions to the eta function 

r 



»J 2V^(^)cosh2 6i V2' + +•''"" / 
sinh6'cosh""^ 6177(5; n+n;/7,„ , B) 



1 r(i + §)„,^ 



^ sinh^ 0cosh^-' ^^ ( 1 1 + 1; 1; _ sinh" ) C ( ^; n+n;/7™, B 



r(^) 



s 3 



2' 2' 



--cosli^-i07,(s;n+n;/7„7;B) 



cosh'*"^ 6* 



sinh6ir(l + f) I s 3 ,2 

_ , ^,0 — Fi[ -;! + s nh^ 

0Fr ^ 2 2' 2' 



C(-;n+n;/7™7,B" 



From here, with the help of H3.o|l and H3.e|l . it is easy to find the residue of 
r]{s; /, P,Bg) at s = to — 1, needed for the evaluation of a^(/, P, Be). We find 



Res ri{m-l-J,P,Be) 



V^r(f) 12 

(to - 1) sinh 9 cosh""^ e 2F1 



— 2-ao(n+n;/7,„,B2 

cosh 9 



1 TO + 1 3 . T 2 /, \ 

-, :-:-smh 6* x 

2' 2 '2' 



[ao(n+n;/7„7, B^) _ tanh0ao(n+n;/7„, B^)] } 
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The leading heat kernel coefficient ao{G,B^) is of course known for an arbitrary 
endomorphism G; it is 



ao(G,B2) = (47r)- 



dy Try (G). 



N 



In order to obtain the invariant form given in Lemma 13.31 we evaluate Il-i-n^ in the 
form 

n+n!^ = - cosh0(cosh0 + 7sinh0 — 77m)- 
Adding up all pieces this shows 

Res ri{m - 1; /, P, Be) = ^^,„^. (47r)-"'/^ / dy Try {hm ■ 



N 



itanhe- i(m- l)sinh6lcosh™"2 6' 2^1 



- , — i — ; - ; — sinh^ 9 
2 2 '2' 



-f-{m~ 1) sinh 6 cosh" 



2F1 



7: 1 — ; 7: ; ^ sinh^ 6 
2' 2 '2' 



-/7 



1 



1 



(m - 1) sinh^ 9 cosh""^ 9 2F1 



2cosh6' 2 

From this result we can read off di{9, m), i = 1, 4; we find 
TO — 1 



7: 1 — ^ — ; 7: ; ~ sinh^ 
2' 2 '2' 



di(6l,m) 

d^{e,m) 
c?4(0, m) 



■ sinh e cosh"" ^61 2^1 



— , — i — ; — ; — sinh^ 9 
2' 2 '2' 



I 



TO — 1 



2 cosh ^ 



■ tanh 



■sinli ^cosh" 



■2F1 



TO — 1 



■sinh6'cosh™-^6' 2F1 



I TO+ f 3 . , 2/, 1 



— , — i — ; — ; — sinh^ 9 
2 2 '2 



The result for di{9,m) can be seen to agree with the result on the ball, eq. p.gj 
by using the transformation formula (|25|. eq. 9.131.1) 



(3.p) 



,Pi(a, /3; 7; ^) = (1 - z)^-"-^ 2^1(7 - a, 7 - /3; 7; 2)- 



In order to show that the results for d2{9,m) coming from the ball and cylinder 
agree, we need to show that 

cosh^^ 2-Fi (^1,1- y;i;-sinh2 
(3.q) = 1 + (m - 1) sinh^ 9 cosh™"^ 9 2^1 (l, |; - sinh" ( 



To see this, we first apply the above transformation formula, eq. p.p| |, and then 
the Gauss recursion formula ([23, eq. 9.137.12) 

7 2Pi(a,/3;7;^) -7 2Pi(a + l,/3;7;z) + /32: 2Fi{a + 1, f3 + + 1; z) = 



with a = -1/2, /? = (m - l)/2, 7 = 1/2, and z = - sinli 6*. Thus, all results 
obtained are consistent and we have determined the full ai and 02 coefficient for 
chiral bag boundary conditions. 
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4. Concluding remarks 

For the case of operators P of Dirac type subject to local boundary conditions of 
chiral bag type as in eq. (1-b), we have studied the asymptotic expansion as t ^ 0^ 
of the smeared L^-trace of the associated heat semigroup, i.e. 

Tr^.(/e-*^^) ^ f;i("-™)/2a„(/,p2^S,). 

n=0 

On using functorial methods, special case calculations and the relation between 77- 
and ^-invariants, we have succeeded in evaluating the full boundary contribution 
to the ai and 02 coefficients, the functional form of which is given by eqs. H2.a|) 
and (2.b). Our contributions are of technical but non-trivial nature, because both 
functorial methods and the theory of the 77-invariant require a lot of work to obtain 
the desired 02 coefficient. It now appears possible that, by exploiting the methods 
described in our paper, further heat-kernel coefficients will be obtained, if they are 
needed in physical or mathematical applications. In turn, a better understanding 
of the spectral functions of modern mathematical physics (23 will also be gained. 
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